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Abstract 

A theoretical analysis of solutions of renormalisation group equations in the MSSM 
corresponding to the quasi-fixed point conditions shows that the mass of the lightest 
Higgs boson in this case does not exceed 94 it 5 GeV. It means that a substan- 
tial part of the parameter space of the MSSM is practically excluded by existing 
experimental data from LEPII. In the NMSSM the upper bound on the lightest 
Higgs boson mass reaches its maximum in the strong Yukawa coupling regime, when 
Yukawa constants are considerably larger the gauge ones on the Grand Unification 
scale. In this paper a particle spectrum in a simple modification of NMSSM which 
leads to a self-consistent solution in the considered region of the parameter space is 
studied. This model allows one to get rrih ~ 125 GeV even for comparatively low 
values of tan/3 > 1.9. For an analysis of the Higgs boson spectrum and neutralino 
spectrum a method for diagonalisation of mass matrices proposed formerly is used. 
The mass of the lightest Higgs boson in this model does not exceed 130.5 it 3.5 GeV. 
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1 Introduction 



The search for the Higgs boson remains one of the top priorities for existing accelerators 
as well as for those still at the design stage. This is because this boson plays a key role 
in the Standard Model which describes all currently available experimental data with a 
high degree of accuracy. As a result of the spontaneous symmetry breaking SU (2) ®U{1) 
the Higgs scalar acquires a nonzero vacuum expectation value without destroying the 
Lorentz invariance, and generates the masses of all fermions and vector bosons. An 
analysis of the experimental data using the Standard Model has shown that there is a 
95% probability that its mass will not exceed 210 GeV At the same time, assuming 
that there are no new fields and interactions and also no Landau pole in the solution 
of the renormalisation group equations for the self-action constant of Higgs fields up to 
the scale Mpi ^ 2.4 x 10^^ GeV, we can show that rrih < 180 GeV 0, 0. In this case, 
physical vacuum is only stable provided that the mass of the Higgs boson is greater than 
135 GeV [0-0. However, it should be noted that this simplified model does not lead to 
unification of the gauge constants and a solution of the hierarchy problem f^. As a 
result, the construction of a realistic theory which combines all the fields and interactions 
is extremely difficult in this case. 

Unification of the gauge constants occurs naturally on the scale Mx ~ 3 x 10^^ GeV 
within the supersymmetric generalisation of the Standard Model, i.e., the Minimal Su- 
persymmetric Standard Model (MSSM) 0. In order that all the fundamental fermions 
acquire mass in the MSSM, not one but two Higgs doublets Hi and H2 must be in- 
troduced in the theory, each acquiring the nonzero vacuum expectation value Vi and V2 
where v"^ = v\ + = (246 GeV)^. The spectrum of the Higgs sector of the MSSM con- 
tains four massive states: two CP-even, one CP-odd, and one charged. An important 
distinguishing feature of the supersymmetric model is the existing of a light Higgs boson 
in the CP-even sector. The upper bound on its mass is determined to a considerable 
extent by the value tan (3 = V2/V1. In the tree-level approximation the mass of the lightest 
Higgs boson in the MSSM does not exceed the mass of the Z-boson {Mz ~ 91.2 GeV): 
< Mz \ cos 2/5 1 0. Allowance for the contribution of loop corrections to the effective 
interaction potential of the Higgs fields from a t-quark and its superpartners significantly 
raises the upper bound on its mass: 



mh<\^Ml cos2 2/? + AiV + ^Xl ■ (1) 



Here A^^^ and A^? are the one-loop |l^ and two-loop [|ll|] corrections, respectively. The 



values of these corrections are proportional to , where is the running mass of t-quark 
which depends logarithmically on the supersymmetry breaking scale Mg and is almost 
independent of the choice of tan/5. In 0, [Q, bounds on the mass of the Higgs boson 
were compared in the Minimal Standard and Supersymmetric models. The upper bound 
on the mass of the light CP-even Higgs boson in the MSSM increases with increasing 
tan /? and for tan /5 ^ 1 in realistic supersymmetric models with M5 < 1000 GeV reaches 
125-128 GeV. 

However, a considerable fraction of the solutions of the system of MSSM renormalisation 
group equations is focused near the infrared quasi-fixed point at tan/5 ~ 1. In the region 
of parameter space of interest to us (tan/? <^ 50) the Yukawa constants of a 6-quark 
(/ife) and a r-lepton {h^) are negligible so that an exact analytic solution can be obtained 
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for the one-loop renormalisation group equations [|T^. For the Yukawa constants of a 
t-quark ht{t) and the gauge constants gi{t) its solution has the following form: 
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bi = 33/5, 62 = 1 
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The variable t is determined by a standard method t = ln(M|-/g^). The boundary 
conditions for the renormalisation group equations are usually set at the grand unifi- 
cation scale Mx {t = 0) where the values of all three Yukawa constants are the same: 
ai{0) = 0:2(0) = 03(0) = (5(0). On the electroweak scale where hf{0) ^ 1 the second term 
in the denominator of the expression describing the evolution of Yt{t) is much smaller 
than unity and all the solutions are concentrated in a narrow interval near the quasi-fixed 
point YQFp(t) = E{t)/6F(t) |jl3[. In other words in the low-energy range the dependence 
of Yt{t) on the initial conditions on the scale Mx disappears. In addition to the Yukawa 
constant of the t-quark, the corresponding trilinear interaction constant of the scalar 
fields At and the combination of the scalar masses 2Jt^ = tuq + mfj + m| also cease to 
depend on At{0) and DJlf{0) as Yt{0) increases. Then on the electroweak scale near the 
infrared quasi-fixed point At{t) and Tl^{t) are only expressed in terms of the gaugino 
mass on the Grand Unification scale. Formally this type of solution can be obtained if 
1^(0) is made to go to infinity. Deviations from this solution are determined by ratio 
l/6F{t)Yt{0) which is of the order of l/10hf{0) on the electroweak scale. 

The properties of the solutions of the system of MSSM renormalisation group equations 
and also the particle spectrum near the infrared quasi-fixed point for tan (3 ~ 1 have 
been studied by many authors [jl4|, |T^. Recent investigations [jl5|-|T^ have shown that 
for solutions Yt{t) corresponding to the quasi- fixed point regime the value of tan/? is 
between 1.3 and 1.8. These comparatively low values of tan/3 yield significantly more 
stringent bounds on the mass of the lightest Higgs boson. The weak dependence of the soft 
supersymmetry breaking parameters At{t) and on the boundary conditions near 

the quasi-fixed point means that the upper bound on its mass can be calculated fairly 



accurately. A theoretical analysis made in [|T5|, [jl6| showed that nih does not exceed 



94 ± 5 GeV. This bound is 25-30 GeV below the absolute upper bound in the Minimal 
Supersymmetric Model. Since the lower bound on the mass of the Higgs boson from 
LEP II data is 113 GeV |jl|, which for the spectrum of heavy supersymmetric particles is 
the same as the corresponding bound on the mass of the Higgs boson in the Standard 
Model, a considerable fraction of the solutions which come out to a quasi-fixed point 
in the MSSM, are almost eliminated by existing experimental data. This provides the 
stimulus for theoretical analyses of the Higgs sector in more complex supersymmetric 
models. 
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The simplest expansion of the MSSM which can conserve the unification of the gauge 
constants and raise the upper bound on the mass of the hghtest Higgs boson is the 



Nonminimal Supersymmetric Standard Model (NMSSM) [jT8|-pO|. In addition to the 
doublets Hi and H2-, the Higgs sector of this model contains the additional singlet su- 
perfield Y relative to the gauge SU{2) (8> U{1) interactions. The most attractive region 
of the NMSSM parameter space from the point of view of theoretical analysis is that 
corresponding to the limit of strong Yukawa coupling when the Yukawa constants on the 
Grand Unification scale Mx are substantially larger than the gauge constant gcvT- This 
is the region where the upper bound on the mass of the lightest Higgs boson reaches its 
maximum, which is several GeVs larger than the corresponding absolute bound in the 
MSSM. In addition, in this particular case it is possible to select the interaction constants 
so as to achieve the unification of the Yukawa constants of a 6-quark and a r-lepton on 
the scale Mx |^ which usually occurs in GUTs [ p3|| . 



However, the mass of the lightest Higgs boson in the NMSSM differs substantially from 



its upper bound pj] . In this connection, the present paper examines a very simple model 
in which rrih reaches its upper theoretical bound for a specific choice of fundamental pa- 
rameters. This bound model is obtained by modifying the Nonminimal Supersymmetric 
Model and yields a self-consistent solution in the strong Yukawa coupling regime where, 
even for comparatively low values of tan/3 < 1.9, the mass of the lightest Higgs boson 
in the modified NMSSM may reach 125-127 GeV. Although the parameter space of this 
model is enlarged considerably, the theory does not lose its predictive capacity. The 
proposed model is used to study characteristics of the spectrum of superpartners of ob- 
servable particles and Higgs bosons. The mass of the lightest Higgs boson in this model 
does not exceed 130.5 ± 3.5 GeV. 

This bound on the mass of the lightest Higgs boson is not the absolute upper bound 



in supersymmetric models. For instance, it was shown in [^] that by introduction four 
or five additional 5-1-5 multiplets of matter, the upper bound on mh in the NMSSM is 
increased up to 155 GeV. Recently the upper bound on the mass of the lightest Higgs 
boson has been actively discussed using more complex expansion of MSSM |2^, p8[] . 
In particular, in addition to the singlet it is also possible to introduce several SU{2) 
triplets into the Higgs sector of supersymmetric models. Their appearance destroys the 
unification of the gauge constants at high energies. In order to reconstruct this, in 
addition to triplets we also need to add several multiplets of matter which carry colour 
charge in the SU{3) group but do not participate in SU{2) U{1) interactions, for 
example four 3-1-3. A numerical analysis made in [^] shows that unification of the 
gauge constants then occurs on the scale Mx ~ 10^"^ GeV and the mass of the lightest 
Higgs boson does not exceed 190 GeV. The existence of a fourth generation of particles 
in the MSSM ||2^, which is extremely problematical from the point of view of the known 



experimental data, also leads to an appreciable increase in the upper bound on m/j. 
Consequently, an increase in the upper bound on the mass of the lightest Higgs boson 
in the supersymmetric models is usually accompanied by a substantial increase in the 
number of particles in the models which may be counted as a serious disadvantage of 



this type of model. In the present study, unlike those noted above |2^, we examine 
the dependence of rrih and the particle spectrum on the fundamental parameters of the 
modified NMSSM in the strong Yukawa coupling regime. 
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2 NMSSM parameters and their renormalisation in 
the strong Yukawa coupling regime 



By definition the superpotential of tlie Nonminimal Supersymmetric Model is invariant 
witfi respect to tlie discrete transformations y'^ = e^'^^^^ya of the Z3 group |19[ which 



means that we can avoid the problem of the /i-term in supergravity models. Zs-symmetry 
usually occurs in string models in which all the fields of the observable sector remain 
massless in the exact supersymmetric limit. In addition to observable superfields ya, 
supergravity theories also contain a hidden sector in which local supersymmetry is broken. 
In modern supergravity theories this sector includes singlet dilaton 5* and moduli fields 
with respect to gauge interactions. These fields always appear in four-dimensional theory 
and they occur as a result of the compactification of additional dimensions. The vacuum- 
averaged dilaton and moduli fields determine the values of the gauge constants on the 
Grand Unification scale and also the dimensions and shape of compacted space. The 
superpotential in supergravity models is usually represented as an expansion in terms of 
superfields of the observable sector p9[|: 



W = Wo{S, Tm) + ^fj^a/siS, Tm)yayp + ^Kp-yiS, T^)yay0y'r, (3) 

where Wo{S,Tm) is the superpotential of the hidden sector. In expression (|^) summation 
is performed over the recurrent greek subscripts. The requirements for conservation of im- 
parity pO| and gauge invariance have the result that the single parameter /i is retained in 



the MSSM which corresponds to the term ^{Hi, H2) in the superpotential @). However, 
the expansion assumes that this fundamental parameter should be of the order of Mpi 
since this scale is the only dimensional parameter characterising the hidden (gravity) 
sector of the theory. In this case, however, the Higgs bosons Hi and H2 acquire an 
enormous mass Tn\^ ~ /^^ ~ "^pi ^"^^ breaking of SU{2) ® U{1) symmetry occurs. 
In the NMSSM the term /j, in the superpotential (^) is not invariant with respect to 
discrete transformations of the Z3 group and for this reason should be eliminated from 
the analysis (/x = 0). As a result of the multiplicative nature of the renormalisation 
of this parameter, the term fi{q) remains zero on any scale q < Mx -7- Mpi. However, 
the absence of mixing of the Higgs doublets on electroweak scale has the result that Hi 
acquires no vacuum expectation value as a result of the spontaneous symmetry breaking 
and d-tjpe quarks and charged leptons remain massless. In order to ensure that all 
quarks and charged leptons acquire nonzero masses, an additional singlet superfield Y 
with respect to gauge SU{2) ® U{1) transformations is introduced in the NMSSM. The 



superpotential of the Higgs sector of the Nonminimal Supersymmetric Model |lq]-[pO| 
has the following form: 

Wh = XY{HiH2) + (4) 

As a result of the spontaneous breaking of SU{2) ® U{1) symmetry, the field Y acquires 
a vacuum expectation value {(Y) = y/V^) and the effective /i-term (/i = \y/^/2) is 
generated. 

In addition to the Yukawa constants A and x, and also the Standard Model constants, 
the Nonminimal Supersymmetric Model contains a large number of unknown parame- 
ters. These are the so-called soft supersymmetry breaking parameters which are required 
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to obtain an acceptable spectrum of superpartners of observable particles form the phe- 
nomenological point of view. The hypothesis on the universal nature of these constants 
on the Grand Unification scale allows us to reduce their number in the NMSSM to three: 
the mass of all the scalar particles mo, the gaugino mass M1/2, and the trilinear inter- 
action constant of the scalar fields A. In order to avoid strong violation of CP-parity 
and also spontaneous breaking of gauge symmetry at high energies (Mpi ^ ^ m^) as 
a result of which the scalar superpartners of leptons and quarks would require nonzero 
vacuum expectation values, the complex phases of the soft supersymmetry breaking pa- 
rameters are assumed to be zero and only positive values of ml are considered. Naturally 
universal supersymmetry breaking parameters appear in the minimal supergravity model 
3T| and also in various string models [^, ||32|. In the low-energy region the hypothesis 



of universal fundamental parameters allows to avoid the appearance of neutral currents 
with flavour changes and can simplify the analysis of the particle spectrum as far as 
possible. The fundamental parameters thus determined on the Grand Unification scale 
should be considered as boundary conditions for the system of renormalisation group 
equations which describes the evolution of these constants as far as the electroweak scale 
or the supersymmetry breaking scale. The complete system of the renormalisation group 
equations of the Nonminimal Supersymmetric Model can be found in |Q, |Q. These 



experimental data impose various constraints on the NMSSM parameter space which 
were analysed in |l36 |. 



The introduction of the neutral field Y in the NMSSM potential leads to the appearance of 
a corresponding F-term in the interaction potential of the Higgs fields. As a consequence, 
the upper bound on the mass of the lightest Higgs boson is increased: 



mh < xI'y^^ sin^ 2(3 + M| cos2 2(3 + A[^^ + A^? . (5) 



The relationship was obtained in the tree-level approximation (An = 0) in . 
However, loop corrections to the effective interaction potential of the Higgs fields from 
the t-quark and its superpartners play a very significant role. In terms of absolute 
value their contribution to the upper bound on the mass of the Higgs boson remains 
approximately the same as in the Minimal SupersjTiimetric Model. When calculating 
the corrections A^^^ and A^^-* we need to replace the parameter jj by Xy/y/2. Studies of 
the Higgs sector in the Nonminimal Supersymmetric model and the one-loop corrections 
to it were reported in [^], |3^-|3^. In ^ the upper bound on the mass of the 



lightest Higgs boson in the NMSSM was compared with the corresponding bounds on 
in the Minimal Standard and Supersymmetric Models. The possibility of a spontaneous 
loss of CP-parity in the Higgs sector of the NMSSM was studied in |]39[, . 



It follows from condition @ that the upper bound on rrih increases as A increases. More- 
over, it only differs substantially from the corresponding bound in the MSSM in the range 
of small tan (3. For high values (tan (3^1) the value of sin 2(3 tends to zero and the upper 
bounds on the mass of the lightest Higgs boson in the MSSM and NMSSM are almost 
the same. The case of small tan (3 is only achieved for fairly high values of the Yukawa 
constant of a t-quark ht on the electroweak scale (/it (to) ^ 1 where to = ln(M|-/m^)), and 
tan (3 decreases with increasing ht{to). However, an analysis of the renormalisation group 
equations in the NMSSM shows that an increase of the Yukawa constants on the elec- 
troweak scale is accompanied by an increase of ht{0) and A(0) on the Grand Unification 
scale. It thus becomes obvious that the upper bound on the mass of the lightest Higgs 
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boson in the Nonminimal Supersymmetric model reaches its maximum on the strong 
Yukawa couphng hmit, i.e., when ht{0) ^ gi{0) and A(0) ^ fi'j(O). 



In our previous two studies ||21|, we analysed the renormalisation of the NMSSM 
parameters in the strong Yukawa coupling regime. We showed |2T[| that as the values of 



the Yukawa constants on the scale Mx increase, the solutions of the renormalisation group 
equations on the electroweak scale are pulled towards a quasi-fixed (Hill) line (x = 0) or 
surface (x 7^ 0) in Yukawa constant space, which limit the range of permissible values of 
ht, A and x. Outside this range in the solutions of the renormalisation group equations 
for Yi{t), where Yt{t) = h1/{A-nf, Yxit) = X^iAnf, and Y^{t) = xV(47r)2, a Landau 
pole appears below the Grand Unification scale and perturbation theory can not be 
applied when ~ M^-. Along the Hill line or surface the values of Yi{t) are distributed 
nonuniformly. As Yi{0) increases , the region in which the solutions of the renormalisation 
group equations are concentrated on the electroweak scale in the strong Yukawa coupling 
regime becomes narrower and in the limit Yi{0) 00 all the solutions are focused near 
the quasi-fixed points. These points are formed as a result of intersection of the Hill line 
or surface with the infrared fixed (invariant) line. This line connects the stable fixed 



point in the strong Yukawa coupling regime |^ with the infrared stable fixed point of 
the system of NMSSM renormalisation group equations The invariant lines their 

properties in the Minimal Standard and Supersymmetric Models were studied in detail 

m. 



m 



As with increasing Yt{0) the Yukawa constants approach the quasi-fixed points, corre- 
sponding solutions for the trilinear constants Ai{t) and the combinations of scalar particle 
masses 97lf(t) 

ml = ml + ml + ml, 
ml = 3m? 



y 



cease to depend on their initial values on the scale Mx- If the evolution of the gauge 
and Yukawa constants is known, the rest of the renormalisation group equations of the 
NMSSM can be considered as a system of linear equations for the soft symmetry breaking 
parameters. In order to solve this system of equations we first need to integrate the 
equations for the gaugino masses and for the trilinear interaction constants of the scalar 
fields Ai{t) and then use the results to calculate DJlj(t). Since the system of differential 
equations for Ai(t) and 3Jl^(t) is hnear, under universal boundary conditions we can obtain 
the dependence of the soft supersymmetry breaking parameters on the electroweak scale 



on A, Ml/2, and ml [^], |46 



A{t) =ei{t)A + fi{t)Mi/2, 
(t) = aiit)ml + hi{t)Ml/2 + Ci{t)AMi,2 + di{t)A^. (6) 

The functions ej(t), fi{t), ai{t), bi{t), Cj(t), and di{t) remain unknown since no analytic 
solution of the complete system of NMSSM renormalisation group equations exists. It 
was shown in that as the quasi-fixed points are approached, the values of the func- 
tions Ciito), ai(to), Ci(to), and di(to) tend to zero whereas for Kj(0) 00 all Ai(t) are 
proportional to M1/2 and all Tlj(t) oc M'^^^- the weak dependence of Ai{t) and Tlf(t) 
in the strong Yukawa coupling regime on the initial conditions has the result that the 
solutions of the renormalisation group equations for trilinear interaction constants and 
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combinations of scalar particle masses and also the solutions for Yi{t) are focused on 
the electroweak scale near the quasi-fixed points. In general under nonuniversal bound- 
ary conditions the solutions for Ai{t) and njt^(t) are grouped near certain lines (x = 0) 
or planes (x 7^ 0) in the soft supersymmetry breaking parameter space. These lines 
and planes are almost perpendicular to the axes At and 9Jt^ whereas the planes in the 
spaces {Af.Ax, A^) and (9^^, 9Jt|, 3Jt^) are also almost parallel to the axes A:^^ and 
Along these lines and planes as l^j(O) increases, the trilinear interaction constants and 
combinations of scalar particle masses go to quasi-fixed points. 



3 Choice of model 

The soft supersymmetry breaking parameters play a key role in an analysis of the particle 
spectrum in modern supersymmetric models. They destroy the Bose-Fermi degeneracy of 
the spectrum in supersymmetric theories so that the superpartners of observable particles 
are substantially heavier than quarks and leptons. However, it should be noted that a 
study of the particle spectrum in the NMSSM is considerably more complex than a study 
of this spectrum in the MSSM for tan/3 ~ 1 since two new Yukawa constants A and 
X appear in the nonminimal supersymmetric model for which the boundary conditions 
are unknown. In turn, the renormalisation of the trilinear interaction constants and the 
scalar particle masses, i.e, the values of the functions ej(to), /j(^o), ai(to), hito), q(^o), 
and di(to), where to = 2\n{Mx/Mf"^^) depends on the choice of ht{0), A(0), and x(0) on 
the grand unification scale. 

The most interesting from the point of view of a theoretical analysis is a study of the 
spectrum of heavy supersymmetric particles when the scale of the supersymmetry break- 
ing M| ^ M|. This is primarily because in this limit the contribution of new particles to 
the electroweak observable ones is negligible (see, for example ^^l)- has been noted, 
the Standard Model highly accurately describes all the existing experimental data. Addi- 
tional Higgs fields and superpartners of observable particles interacting with vector 
and Z bosons make a nonzero contribution to the electroweak observables. However, for 
M| ^ M| their contribution is suppressed in a power fashion as {Mz/MsY ■, where any 
increase in the scale of the supersymmetry breaking leads to convergence of the theo- 
retical predictions for the strong interaction constant as{Mz) which may be obtained 



assuming unification of the gauge constants with the results of an analysis of the 



experimental data In addition, it should be noted that the mass of the hghtest Higgs 
boson which is one of the central objects of investigation in any supersymmetric model 
reaches its highest value for Ms ~ 1 — 3 TeV. 

Unfortunately, in the strong Yukawa coupling regime in the NMSSM with a minimal 
set of fundamental parameters it is impossible to obtain a self-consistent solution which 
on the one hand would lead to a spectrum with heavy superparticles and on the other 
could give a mass of the lightest Higgs boson greater than that in the MSSM. When 
calculating the particle spectrum, the fundamental parameters A, rriQ, and M1/2 on the 
scale Mx should be selected so that the derivatives of the interaction potential of the 
scalar potential of the scalar fields V{Hi, H2, Y) with respect to the vacuum expectation 
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values Vi, V2, and y would be zero at the minimum: 



dV{vi,V2,y) ^ Q dV{vi,V2,y) ^ ^ dV{vi,V2,y) ^ ^ 
dvi ' 9f2 ' dy 

Since the trilinear interaction constants and the scalar particle masses in the strong 
Yukawa coupling regime are almost independent of A, Eqs.(|^ link the vacuum expecta- 
tion value of the neutral scalar field (Y), and the parameters rriQ and M1/2. The value 
of tan f3 is determined using the Yukawa constant of a t-quark on the electroweak scale 
(see below). Then a spectrum of heavy supersymmetric particles is only achieved when 
A/x ^ 1. However, in this region of parameter space the value of becomes negative 
and the physical vacuum is unstable which can be attributed to the strong mixing of 
the CP-even components of the neutral field Y and the superposition of Higgs doublets 
h = Hi cos (3 + H2 sin [3. 

Studies of the particle spectrum in the nonminimal supersymmetric model 
3^,[^,[^,[Q have shown that a self-consistent nontrivial solution of the system of 



nonlinear algebraic Eqs.(0) for \ (Y) \ < 10 TeV which determines the position of the min- 
imum of the interaction potential of the scalar fields, only exists for A^(to), x^(^o) ^ 0.1. 
In this case a strict correlation exists between the fundamental parameters of the NMSSM. 
In particular, in order to ensure that spontaneous breaking of SU{2) ® ^(1) symmetry 
occurs and the field Y has a nonzero vacuum expectation value of the field, the condition 
\A^/my\ > 3 must be satisfied. However, the following inequalities must also be satisfied: 

Af <3imj + ml^+mlj, 
Al<3{ml + ml^+mlJ, 
Al<3iml + ml^+mlJ. 

Otherwise, the superpartners of leptons and quarks acquire vacuum expectation val- 
ues All these constraints have the result that the ratio |y4/mo| varies between 3 
and 4. In |5^, the particle spectrum in the NMSSM is analysed separately for 
tan/? = mt{rrit) / rrtbimt) and under nonuniversal boundary conditions. 

The limit ^ A^, in the nonminimal supersymmetric model corresponds to the 



MSSM 1^31 . For x = the Lagrangian of the Higgs sector of the NMSSM is invariant 
with respect to the global SU{2) ® U{1) ® U{1) transformations. As a result of the spon- 
taneous symmetry breaking, only the U{1) symmetry corresponding to electromagnetic 
interaction remains unbroken, which leads to four massless degrees of freedom. Two of 
these are eaten by a charged boson and one by a Z boson. Ultimately, the spectrum 
of the nonminimal supersymmetric model for x = contains one physical massless state 
which corresponds to the CP-odd component of the field Y . For low values x^ ^ -^^^t 
the mass of the lightest CP-odd boson is nonzero and is proportional to the self-action 
constant of the neutral superfield Y . If the Yukawa constants A, x ~ 10"'^ — 10~^, for a 
certain choice of fundamental parameters the mass of the lightest CP-even Higgs boson 
may be only a few GeV [^,^^,[^. The main contribution to its wave function is made 



by the neutral scalar field Y which makes it very difficult to search for this on existing 
accelerators and those at the design stage since the interaction constants of this type of 
Higgs boson with gauge bosons and fermions are small. In this limiting case, the lightest 
stable supersymmetric particle having R-parity of —1 is usually the superpartner of the 
neutral scalar field F BsH, i6|. 
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However, unlike the minimal supersymmetric model, the discrete Z3 symmetry which 
can avoid problems of the /x-term in the NMSSM has the result that three degenerate 
vacuums appear in the theory because of the breaking of gauge symmetry. Immediately 
after a phase transition on the electroweak scale the is filled equally with three degenerate 
phases. The entire space is then divided into separate regions in each of which a partic- 
ular case is achieved. The regions are separated by domain walls with the surface energy 
density a ^ . Data from cosmological observations eliminate the existence of domain 
walls. The domain structure of vacuum in the NMSSM is destroyed if the vacuum degen- 
eracy caused by symmetry disappears. It was shown in that breaking of 
symmetry by introducing into the NMSSM Lagrangian nonrenormalisable operators of 
dimension d = 5 which do not break the SU{2) ® U{1) symmetry can be used to obtain 
splitting of initially degenerate vacuums such that the domain walls disappear before 
the beginning of the nucleosynthesis are (T ~ 1 MeV). Although operators of dimension 
d = 5 are suppressed with respect to Mpi in supergravity models, their introduction leads 
to quadratic divergences in the two-loop approximation, i.e., to the problem of hierar- 
chies. Consequently, linear and quadratic terms with respect to superfields are generated 
in the superpotential of this theory and the vacuum expectation value of the neutral 
scalar Y is of the order of 10^^ GeV. 

In order to avoid a vacuum domain structure and obtain a self-consistent solution in 
the strong Yukawa coupling regime, we need to modify the nonminimal supersymmetric 
model. The NMSSM can be modified most simply by introducing additional terms in 
the superpotential of the Higgs sector: fi{HiH2) and fi'Y"^ which are not forbidden by 
gauge SU{2) U{1) and R symmetries. The additional bilinear terms in the NMSSM 
superpotential destroy the Z3 symmetry and no domain structures appear in this theory 
since no system of degenerate vacuums exists. The introduction of the parameter fi 
ensures that it is possible to obtain a spectrum of heavy supersymmetric particles in 
the strong Yukawa coupling regime in the modified model and for a certain choice of fi' 
the mass of the lightest Higgs boson reaches its upper bound. In this case the mass of 
the lightest Higgs boson has its highest value x = since as x(to) increases, the upper 
bound on its mass is reduced as a result of a decrease in X(to). In the limit x = the 
CP-odd Higgs sector of the modified NMSSM contains no physical massless states since 
in this case, the global symmetry of the Lagrangian is the same as the local symmetry 
which eliminates the Yukawa self-action constant of the neutral field Y from the analysis. 
Assuming that this is zero and neglecting all the Yukawa constants except for A and ht, 
the complete superpotential of the modified NMSSM can be expressed in the following 
form: 

W^NMSSM = KH1H2) + /i'r^ + XY{HiH2) + htiH^QM ■ (8) 

where f/^ is the charge-coupled right superfield of a t-quark and Q is a doublet of left 
superfields of b and t-quarks. 

In supergravity models, bilinear terms with respect to the superfields may be generated 
in the superpotential (§) as a result of the additional term Z{HiH2) + h.c. in the Kahler 
potential [^,|5^ or nonrenormalisable interaction of the fields of the observable and 
hidden sectors. The appearance of nonrenormalisable operators of this type in the su- 
perpotential of supergravity models may be attributed to nonperturbative effects (for 
instance, gaugino condensation) [Q,||59|. In addition to the parameters /x and /i', this 



model also sees the appearance of the corresponding bilinear interaction constants of the 
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scalar fields B and B' which for a minimal choice of fundamental parameters should be 
assumed to the equal on the grand unification scale. Thus, the nonminimal supersym- 
metric model may include seven fundamental parameters in addition to the constants of 
the Standard Model: 

A, fx, n', A, B, mo, M1/2 . 



4 Constraints on the parameter space of the modi- 
fied NMSSM 

Despite a substantial expansion of the parameter space, the theory does not lose its 
predictive capacity. An analysis of the behaviour of the solutions of the NMSSM renor- 
malisation group equations in the strong Yukawa coupling limit for x = showed that 
for Yi{0) —>■ 00 all the solutions are concentrated near the quasi-fixed point: 

pf^'ito) = 0.803 , pfito) = 1.77, pjf (to) = 6.09 

* (9) 
pr(to) = 0.224 , pZ^'Ho) = -0.42 , pg,f (to) = -2.28 , 

where p,(t) = ^ p.(t) = ^ p^.(t) = ^ and p^.(t) = 
a3(t) a3(t) Ml/2 ^ 

Thus, at the first stage of the analysis we fixed the initial values of the Yukawa 

constants A^(0) = hf{0) = 10 corresponding to the quasi- fixed point regime @ 

of the renormalisation group equations and also the supersymmetry breaking scale 

M3(1000 GeV) = 1000 GeV which determines the mass scale of all the supersymmet- 

ric particles. 

Existing FNAL experimental data from measurements of the mass of a t-quark can 
uniquely relate tan/5 to the Yukawa constant ht of a t-quark. The running mass of a 
t-quark generated when the SU{2) ® U{1) symmetry is broken is directly proportional 
to /it(to): 

mt{Mr'') = ^'^^PK sin f3. (10) 
v2 



However, the value of mt{Mf"^^) calculated in the MS scheme is equal to 

'mt{Mf°^'^) = 165 ± 5 GeV. The inaccuracy in determining the running mass of a t- 
quark is primarily attributable to the experimental error with which its pole mass is 
measures {Mf"^^ = 174.3 ±5.1 GeV [^). For each fixed set of boundary conditions ht{0) 
and A(0), using renormalisation group equations we can calculate the Yukawa constant 
of a t-quark on the electroweak scale and then, substituting the value obtained /it(to) 
into formula (|TD|), we can determine the value of tan/?. In the infrared quasi-fixed point 
regime we obtain tan/3 ~ 1.88 for mt{Mf°^^) = 165 GeV (Tables 1 and 2). 

An additional constant which fixes M3 can be used to determine one of the supersymmetry 
breaking parameters Mi/2- The values of all the other dimensional parameters p, p', A, B, 
and niQ should be selected so that spontaneous breaking of gauge SU (2) (S)t/(1) symmetry 



11 



occurs on the electroweak scale. The complete interaction potential of the Higgs fields in 
the modified NMSSM can be expressed as the sum: 



V{H^, H2, Y) = ^Ii\H,\' + fii\H2\' + fiy\Y\' + [fii{H,H2) + fiiV + XAyXiH.H^) + 
+ Xfi'Y*{HiH2) + XfiY{\Hi\^ + \H2\^) + /i.e.] + A'Kifiifa)!' + X^Y^Hil^ + \H2\^) + 

+ ^{\H2\'- \H,\Y + -^{Ht(TH, + HtfTH^f + ^V{H,, H^, Y) , (11) 

o o 

where AV{Hi, H2,Y) are the one-loop corrections to the effective interaction potential; 
g and g' are the constants of the gauge SU{2) and U{1) interactions (^^i = ^^/5/3g'). 
The constants in the interaction potential ( |Tl|) are related to the soft supersymmetry 
breaking parameters as follows: 

2 2,2 2 2,2 2 2|/2 

^1 = + /i , = m2 + /i , fiy = my + , 

where 

ml{Mx) = ml{Mx) = ml{Mx) = ml , 
BiMx) = B'iMx) = Bo. 

The position of the physical minimum of the interaction potential of the Higgs fields ([TT| ) 
is determined by Eqs.(0). Since the vacuum expectation value v and tan (3 are known, 
the system of Eqs.(0) can be used to find fi and Bq. Then, it is convenient to introduce 

Hes = + instead of After various transformations we obtain: 



2 _ ml - ml tan^ /3 + Azjfics) _ 1 , ,2 
tan2/5-l 2^^^^ 



+ + 2fils + + A/3(m) sin 2(3 = -2 ("e/x + ) (12) 



y (ml + fj.'^ + B'fi'^ = ^v^X, - A,(/ieffj 



where Aj corresponds to the contribution of the one-loop corrections: 

- 2 dAV dAV _ dAV 

^ ~ vHan2(3 d(3 ^ ' ' ~ ~df ' 

I { dV 1 dAV 

Az = —1^ 2— cos 2/? - -— — sin 2(3 
cos^p I ov^ op 



and 



Xi = (2/ieff + (/i' + Ax) sin 2p) , X2 = (/i' + Aa) cos 2/3. 

When calculating the one-loop corrections we shall only take into account the contri- 
bution from loops containing a t-quark and its superpartners since their contribution is 
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dominant. In supersymmetric theories each fermion state with a specific chirahty has 
a scalar superpartner. Thus, a t-quark incorporating left and right chiral components 
has two scalar superpartners, right t/j and left t^, which become mixed as a result of the 
spontaneous breaking of SU{2) ® U{1) symmetry, and this results in the formation of 
two charged scalar particles having masses m^^ and mi^: 



= \i^Q + ^l + 2m? ± ^{ml-mlf + AmiX^) , (13) 

where Xt = At + /icfi/tan/?. Since m| and m| should be positive we have 

Xf < — 2^^Q ~^ ^t)i^u + ^t)- Otherwise, the quark fields acquire nonzero vacuum 
rrit 

expectation values and the gauge SU{2) ® U{1) symmetry of the initial Lagrangian is 
completely broken, which leads to the appearance of nonzero masses for gluons and pho- 
tons. The contribution of the one-loop corrections from the t-quark and its superpartners 
to the effective interaction potential of the Higgs fields is expressed only in terms of their 
masses: 




(14) 



For this reason all Aj are merely functions of /ies and do not depend on Bq and y. 



Using the first equation of the system ([T^) we can find /ieg. In this case, the sign of /ieg is 
not fixed and must be considered as a free parameter in the theory. Substituting this value 
of /ieff into the two remaining equations of the system (|12D , we can eliminate Bq from the 
number of independent fundamental parameters and calculate the vacuum expectation 
value of the field Y: {¥) = y/V^ and in order to find Bq we need to bear in mind the 
relationships linking the bilinear interaction constants B and B' in the electroweak scale 
with Bq obtained by solving the renormalisation group equations in the modified NMSSM 
(see Appendix): 

B{t) = at)Bo + (T{t)A + cu(t)Mi/2, 

B'{t) = Bo + ai{t)A^ + ^i(t)Mi/2^ , (15) 

where (i't), c"(^); (^lit), uj{t), and LJi(t) are various functions of t, ht{0), and A(0) which 
do not depend on the choice of fundamental soft supersymmetry breaking parameters on 
the grand unification scale, and also on /xq and /Xg. For a fixed sign of fies three different 
solutions of the system of Eqs.([I^) exist. However, only one of these is of interest from 



the physical point of view. It follows from the last equation in ([T^ ) that the vacuum 
expectation value of the field Y is of the order of ?/ ~ Xv'^/Ms and for the case of heavy 
supersymmetric particles y <^ v. The other two solutions give an excessively light CP- 
even Higgs boson which corresponds to fine tuning between the fundamental parameters 
Bq and /i'. 

The parameters fies and Bq thus determined and also the vacuum expectation value y 
depend on the choice of A, mg, and /i'. Thus, at the next stage of the analysis of the 
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modified NMSSM we studied the dependence of tlie particle spectrum on tliese funda- 
mental parameters using Eqs.(P) linking v4j(to) and 9Jt^(to) to A and mg. Similarly we 
investigated the spectrum of superpartners of observable particles and Higgs bosons for 
other values of the Yukawa constants from the vicinity of the infrared quasi- fixed point. 
Although for tan < 2 in the strong Yukawa coupling regime the parameters ht and A can 
be selected so that the Yukawa constants of a 6-quark and a r-lepton would be the same 
on the grand unification scale |^ , , when studying the particle spectrum in the mod- 
ified NMSSM we do not confine ourselves to the case Rbri^) = 1 where Rbr = hb{t) /hr{t). 
This condition arises in minimal schemes of gauge interaction unification |^ and imposes 
very stringent constraints on the parameter space of the model being studied. However, 
since hf, and h^- for tan /3 ~ 1 have small absolute values, they can be generated by means 
of nonrenormalisable operators as a result of the spontaneous symmetry breaking on the 
scale Mx and in this case, the Yukawa constants of a 6-quark and a r-lepton may differ. 



5 Calculations of masses of Higgs bosons and neu- 
tralinos 

We shall first consider the Higgs sector in the modified NMSSM which includes three 
CP-even states, two CP-odd, and one charged Higgs boson. The determinants of the 
mass matrices of the CP-odd and charged Higgs bosons go to zero which corresponds 
to the appearance of two Goldstone bosons which are eaten by massive vector and 
Z bosons during spontaneous breaking of SU{2) (g) U{1) symmetry. The (3 x 3) mass 
matrix of the CP-odd sector is formed by mixing the imaginary parts of the neutral 
components of the Higgs doublets with the imaginary part of the field Y. However, since 
the determinant of this matrix is zero, the problem of finding the eigenvalues reduces to 
solving an ordinary quadratic equation. The calculated masses of the CP-odd states in 
the modified NMSSM are 




= ^\ m\ + <±\lim\- m\f + 4 ( ^(/i' + A) + Ao ) | , (16) 



2 



X2 \ 

2 2, 2io2i^2iA 2 2|/2 d//i 2iA 

= + m2 + + + Aa, = + /i - B ii ^ —v + A3, 

where Aj are the one-loop corrections [l33| , [0^ , |]38 



A more complex situation is encountered in the sector of the CP-even Higgs fields which 
appears as a result of mixing of the real parts of the neutral components of the Higgs 
doublets and the field Y . The determinant of the mass matrix of the CP-even sector ix 
nonzero and thus in order to calculate its eigenvalues we need to solve a cubic equation. 
However, for the case of heavy supersymmetric particles [Ms ^ Mz) in the Higgs field 
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basis 



V 2 V 2 



X2 
X3 



sin /3 Re + ^ cos/5 Re 
V2 V2 



(17) 



1 



V2 



ReY 



this matrix has a hierarchical structure: 








0^ 






El 





H 















Vl2 


Vl3 






V23 


VV^31 


Vs2 


^33 



where 



E! = 0, 



1 



eflf) 



Vn = M| cos' 2/5 + -X'v' sin" 2/5 + A 



11) 



V^i3 = V^3i = At;Xi + Ai3 



V^12 = V, 



21 



^AV -^M|) Sin 4/? + A12, 



1^22 = Mi sin' 2P + -X'v' cos' 2/5 + A^ + A22, 



V22. = V32 = Af X2 + A23 , 
1 



33 



-AV + A33 . 



Here A^ and Ajj are the corrections from loops containing the t-quark and its super- 
partners. The hierarchical structure of the mass matrix means that perturbation from 
quantum mechanics can be used to diagonalise it. The role of the smallness parameters 
in the perturbation theory are played by the ratios M'^/E^ and M'^/E^. This method of 



calculating the masses of Higgs bosons in supersymmetric theories was developed in |2| 
Also discussed there is the simplest method of obtaining a hierarchical mass matrix in the 
Higgs field basis {^7^- A numerical analysis made in this study showed that perturbation 
theory can be used to calculate the masses of Higgs bosons in the modified NMSSM to 
within 1 GeV (~ 1%). 

This method can be used to diagonalise the mass matrix of a neutralino which occurs 
as a result of the mixing of superpartners of gauge bosons W3 and B (or Z and 7) with 
superpartners of neutral Higgs fields. In the basis {B^, W'^, H2, Y), this mass matrix 
has the following form: 

\ 


^ • 

—j=v sm p 
-^v cos p 



/Ml 





-A' 


B' 





M2 


C 


-D' 


-A' 


C 







B' 


-D' 








v 



—=vsmp —=vcosp 
V2 V2 



(19) 



The fragment of the (4 x 4) matrix which includes the first four columns and for rows is 
the same as the mass matrix of a neutralino in the MSSM with A', B', C, and D' given 



15 



by 



A' = Mz cos P sin 9w, 
C = Mz cos P cos 9w, 



Using the unitary transformation U: 



U 



(\ 

1 





\0 






1 



B' = Mz sin (3 sin Ow, 
D' = Mz sm (3 cosOw 






1 

? 



V2 V2 









1/ 



the matrix (^) can be reduced to the form (|T^) M' = UMU^ and then using the ratios 
(M^Z/ieflf)^ and [Mz/fi'^ as the small parameters in the first order of perturbation theory 
for the spectrum of supersymmetric particles we obtain 



m 



XI 



m 



X2 



m 



X3 



m 



X4 



"^X5 



Mi + 
M2 + 







(A' + 


2 (Ml 




^ 2(Mi+/ieff) ' 


(C- 




(C + D'f 
+ 2(M2 + M ' 


2(M2 


- /^cff) 




B'f 


(C - D'f 

+ 2(/ieff-M2) ^ 


2(Aicfr 


-Ml) 



/\2 



{A' + 

2 (Ml + /Xeff) 2(M2 + M 



X^v^sin\(3 + Ti/A) 
2(/icfr - /i') 
(C' + D')' AVsin2(/?-7r/4) 



(20) 



2(/Xeff + A^O 



, AVsin2(/? + 7r/4) AV sin2(/5 - vrM) 
H ^ \- 



2(/i' - /icflf) 



2(/i' + /icff) 



The accuracy with which iriy^- is calculated is slightly lower than that for the CP-even 
Higgs sector. This primarily because the parameters used for the expansion when diago- 
nalising the neutralino mass matrix according to perturbation theory are larger. At this 
point we shall not discuss the spectrum of squarks, sleptons, and charginos in greater 
detail since the analytic expressions for the masses of these particles remain the same as 
in MSSM. In this case, in all the formulas we need to replace /i with /igff. We merely note 
that in the principal approximation the masses of two Dirac charginos and neutralinos 

M2. 



X2 and Xs are the same: m-± 



Atcff and m-± 



6 Results of the numerical analysis 

Results of a numerical analysis of the spectrum of Higgs bosons and superpartners of 
observable particles in the modified NMSSM are given in Figs. 1-5 and Tables 1-3. We 
first need to note that for a fixed sign of yUefr there are two allowed regions of parameter 
space. In one of these the mass of the lightest Higgs boson is greater than in the MSSM 
(see Figs, la and 3a) whereas in the other it is smaller (see Figs, lb and 3b). The mass 
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of the lightest calculated in the first order with respect to perturbation theory has the 
following form: 

ml^Vu-^^- Ml cos^ 2, + \X^.^ s,„^ 2, + A„ - ^^^^^f^ • (21) 

Since the matrix element V12 ~ M|, we neglected its contribution to irih. 
The mass of the lightest CP-even Higgs boson reaches its highest value when 

, _ AMr _ A _ — 13 ]^ Y _ Q /• Figs, la and 3a). Thus, nih is 

^ sin 2/3 Aw sin 2/3' v & j , h 

larger in that region of parameter space where the signs of /i' and fics are opposed. In the 
limit n' ±00 the masses of the CP-even and CP-odd Higgs bosons corresponding to 
the field Y become much larger than the scale of the supersymmetry breaking. In the low 
energy region their contribution to the effective interaction potential of the Higgs fields 
Hi and H2 disappears and the mass of the lightest Higgs boson is the same as in the 
MSSM. For this reason, as can be seen from the graphs plotted in Figs, la, lb and 3a,3b, 
ruh reaches a constant value when fi' — > ±00. The results of the numerical calculations 



plotted in these figures indicate that the two-loop corrections [TT[ play a significant role 
in the calculations of the mass of the lightest Higgs boson. In this particular case, they 
reduce its mass approximately by 10 GeV. Although the one-loop corrections increase 
logarithmically as the scale of the supersymmetry breaking increases, their increase for 
Ms ^ Mz is completely compensated by the log-log asymptotic form of the two-loop 
corrections and rrih remains almost constant. Allowance for the loop corrections has the 
result that for fi^s < the mass of the lightest Higgs boson is greater than that in the 
case /Xeff > 0. This can be attributed to the fact that rrih increases with increasing mixing 
in the superpartner sector of the t-quark (Ir and ii) which is determined by the value 
of Xt = At + fies/ tan/3. Since At < 0, the absolute value of the mixing between in and 
ii is greater for ^^s < 0. It should be noted that the mass of the lightest Higgs boson is 
almost independent of A and mo because of the weak dependence of the squark mass on 
the corresponding fundamental parameters (see Tables 1 and 2). 

Since that part of the parameter space in which fies and /i' have the same sign is almost 
eliminated by the existing experimental data, it is most interesting to study the spectrum 
of Higgs bosons in the region where the mass of the lightest Higgs boson is greater than 
that in the minimal supersymmetric model. In this particular region of parameter space 
the bilinear soft supersymmetry breaking constants and fies have opposite signs and near 
the maximum of ruh the parameter is /i' ~ ~2^cs/ sin 2/3 > \fics\ ~ Ms). For this 
reason the heaviest particle in the modified NMSSM spectrum is the CP-even Higgs boson 
which corresponds to the neutral field Y sice its mass in the principal approximation with 
respect to perturbation theory is M| ^ > /i'^ and is substantially larger than the 
scale of the supersymmetry breaking. It can be seen from Figs. 2a and 4a that the 
mass of the other heavy CP-even Higgs boson {rriH) is almost independent of fi' since 
m| ^ rnjj. However, the spectrum of the CP-odd Higgs sector is determined to a 
considerable by the choice of fundamental parameters. As fi' increases, the mass of the 
CP-odd Y increases and the latter becomes one of the heaviest particles. For values of 
fi' ~ B' the mass of the lightest CP-odd Higgs boson itl\^ ^ is very low (see Figs. 2a, 
4a), which leads to the appearance of a constraint on /x'. Nevertheless, with this choice 
of fundamental parameters this Higgs boson is negligibly involved in electromagnetic and 
weak interactions since the main contribution to its wave function is made by the CP- 
odd component of the field Y. Thus, even when its mass is relatively low, it is extremely 
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difficult to detect this particle experimentally, the heaviest fermion in this model is the 
superpartner of the field Y. Its mass m^^ is proportional to fi' (see (|20|) ). The spectrum 
of remaining neutralinos, charginos, squarks, and sleptons does not depend on the choice 

of fl'. 

Since the dependence of the soft supersymmetry breaking parameters on A disappears in 
the strong Yukawa coupling regime on the electroweak scale, the spectrum of superpart- 
ners of the observable particles and also Hes and B whose numerical values are determined 
by solving the system of Eqs.(^), vary weakly when the trilinear interaction constant 
of the scalar fields varies between — M1/2 and Mi/2- Despite this, the dependence of the 
Higgs boson spectrum on A is conserved. This is mainly because the bilinear interaction 
constant of the neutral scalar fields B' is proportional to A. Using the relations (p3D, we 
obtain 



Ml/2 



where x = A/M1/2. As x increases for < {fies > 0) the bilinear interaction constant 
decreases (increases) in absolute value and conversely. As \B'\ decreases, the masses of 
the CP-even and CP-odd states corresponding to the neutral field Y converge. At the 
same time, an increase in the absolute value of B' leads to a decrease in which 
disappears when B' ~ fi'. The dependence of the Higgs boson spectrum of the parameter 
A for mo = is studied in Figs. 2b and 4b. The parameter fi' in this particular case is 
selected so that the mass of the lightest Higgs boson coincides with the upper bound on 
ruh for A = 0. 

Although in some cases we assumed mo = when analysing the modified NMSSM, this 
limit is unacceptable from the physical point of view since in this case the lightest (and 
consequently stable) supersymmetric particle is the superpartner of the right r-lepton 
which contradicts existing astrophysical observations. However, as mo increases the mass 
of the superpartner of the right r-lepton increases and even for comparatively low values 
of mo/Mi/2 the lightest particle in the spectrum of superpartners of observable particles 
becomes the neutralino. The results of the numerical calculations presented in Tables 
1 and 2 can be used to assess the infiuence of the fundamental constants A, mo, and 
Ml/2 on the superpartner spectrum of the t-quark, gluinos, neutralinos, charginos, and 
Higgs bosons. For each set of parameters listed above we give the values of the upper 
bound on the mass of the lightest Higgs boson, calculated in the one-loop and two-loop 
approximations and also the corresponding fi^s, Bq, y, and fi' for which V13 = 0. It can 
be seen from the data presented in Table 1 that the qualitative pattern of the spectrum 
remains unchanged if the parameters A and mo vary within reasonable limits. It should 
also be noted that as ml increases, the masses of squarks, sleptons, Higgs bosons, and also 
heavy charginos and neutralinos increases whereas the spectrum of the lightest particles 
remains unchanged. The mass of a charged Higgs boson which has not been mentioned 
before is almost independent of A and fi' and numerically similar results are obtained for 
the mass of a charged Higgs boson mH±- 

In the present paper we have made a detailed study of the superpartner and Higgs boson 
spectrum for initial values of the Yukawa constants /i^ (0) = A^(0) = 10 corresponding to 
the scenario of the infrared quasi-fixed point in the NMSSM. The results of the numerical 
calculations presented in Tables 1 and 2 indicate that for mt(Mf°'^) and < 2 TeV the 
mass of the lightest Higgs boson does not exceed 127 GeV. Other data presented in Table 



18 



3 indicate that the distinguishing features of the supersymmetric particle spectrum are 
conserved for h^{0) ^ A^(0) and hf{0) <^ A^(0) as along as the Yukawa constants on the 
grand unification scale are substantially larger than the gauge constants. Nevertheless, 
the upper bound on the mass of the lightest Higgs boson, the value of tan/?, and the 

particle masses calculated for u' = — Ax — — when V13 = vary as a 

sm 2p Xv sm 2p 

function of the choice of h^{0) and A^(0). Nevertheless, as A^(0) decreases from 10 to 2, 

the upper bound on for M3 = 1 TeV drops from 128 to 113 GeV (see Table 3). Thus, 

at the concluding stage of the analysis of the modified NMSSM for each fixed tan (3 we 

selected the Yukawa constant A (to) so that reached its highest value on condition that 

perturbation theory can be applied as far as the grand unification scale. The dependence 

mft(tan P) thus obtained is plotted in Fig. 5 where we also plotted the upper bound nih in 

the MSSM as a function of tan f3. As was to be expected, the two bounds on the mass of 

A^w^ 

the lightest Higgs boson are almost the same for large tan P when the term — ^ — sin^ 2/? 

in Eq.(|^) tends to zero. The curve m/i(tan/5) in the NMSSM reaches its maximum when 
tan/5 ~ 2.5 which corresponds to the strong Yukawa coupling regime. Both bounds on 
the mass of the lightest Higgs boson were obtained for M3 < 2 TeV. By varying the scale 
of supersymmetry breaking we can show that the mass of the lightest Higgs boson in the 
NMSSM does not exceed 130.5 ± 3.5 GeV. The indeterminacy observed in calculations 
of the upper bound on m/j is mainly attributable to the experimental error with which 
the mass of a t-quark is measured. 



7 Conclusions 

In the nonminimal supersymmetric model the mass of the lightest Higgs boson reaches its 
highest value in the strong Yukawa coupling regime when all the solutions of the renor- 
malisation group equations are grouped near the infrared quasi-fixed point. However, in 
this region of the parameter space using the NMSSM with a minimal set of fundamental 
parameters it is not possible to obtain a self-consistent solution which on the one hand 
would give a spectrum with heavy supersymmetric particles and on the other could give 
a mass of the lightest Higgs boson greater than that in the MSSM. In order to find such 
a solution, we need to modify the nonminimal supersymmetric model. In the present 
paper we studied the spectrum of superpartners and Higgs bosons using a very simple 
expansion of he NMSSM which can give a self-consistent solution in the strong Yukawa 
coupling regime. Although the parameter space of this model is expanded substantially, 
the theory does not lose its predictive capacity. 

The mass matrix of the CP-even Higgs sector in the modified NMSSM has a hierarchical 
structure which means that it can be diagonalised using a method of calculating the 
spectrum of Higgs bosons proposed earlier, which is based on the ordinary perturbation 
theory of quantum mechanics. This method can be used to calculate the mass of Higgs 
bosons to within 1 GeV (~ 1%). In this case the mass of the lightest Higgs boson near 
the infrared quasi-fixed point for mt{Mf"^^) = 165 GeV and M3 < 2 TeV does not exceed 
127 GeV. By varying the ratio of the Yukawa constants on the grand unification scale, we 
can show that rrih < 130.5 ±3.5 GeV where the indeterminacy observed when calculating 
the upper bound on m/j is mainly attributable to the experimental error with which the 
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mass of the t-quark is measured. The heaviest particle in the region of the parameter 
space of interest is the CP-even Higgs boson corresponding to the neutral field Y. 

In the present study we used the same method of diagonalising the mass matrices to 
calculate and study neutralino masses. As a result we showed that the heaviest fermion 
in the dominant region of parameter space is Y, the superpartner of the neutral scalar field 
Y. For values of < M^^^ gluinos, squarks, heavy CP-even and CP-odd Higgs bosons 
are substantially heavier than sleptons, lightest charginos, and neutralinos. The only 
exception is one of the CP-odd Higgs bosons whose mass varies substantially depending 
on the choice of parameters of the model. However, even if it is relatively low, for example, 
of the order of Mz, there are certain problems involved in recording it experimentally 
since the main contribution to its wave function is made by the CP-odd component of 
the field Y. 

The upper bound on the mass of the lightest Higgs boson in the nonminimal super- 
symmetric model was also studied on recent publications pSf and The predictions 



obtained in these studies are 5 — 6 GeV higher than the bound given above. The dif- 
ference in the predictions can be attributed to the fact that the authors of and 



used the value of \Xt/Ms\ = where Ms = ^/m^JTn^ to calculate the upper bound on 
m/j since the mass of the lightest Higgs boson reaches its highest value for this value of 
Xf. However, in the strong Yukawa coupling regime in the modified NMSSM the ratio 
\Xt/Ms\ is 1.4 — 1.5. Since the mass of the lightest Higgs boson increases with increasing 
mixing between the t-quark superpartners for < \Xt/Ms\ < VQ, and the ratio \Xt/Ms\ 
is considerably less than y^, the upper bound on ruh in the realistic expansion of the 
NMSSM is more stringent that the absolute bound in the nonminimal supersymmetric 
model. 
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Appendix 



Renormalisation group equations for /i, /i', B, and B' pctrameters in the mod- 
ified NMSSM and their solutions. 

Side by side with the trihnear couphngs and masses of scalar particles the modified 
NMSSM, within which we investigate the particle spectrum, contains the parameters 
/i', B, and B'. An evolution of these constants is described by four renormalisation group 
equations: 

-^--i2Y, + 3Y,-3a,--a, 



d4j!_ 

'dt 
dB 

'dt 

dB' 

~dt 



\ H 5 

- (2Y^B' + 4^1^54 + 4Fa^a4 + 41^-^-1 • 



(22) 



For X = with the minimal set of fundamental parameters -B(O) = -B'(O) = Bq, 
Ai{0) = A, Mj(0) = Ml/2, /^'(O) = /Xq, and /x(0) = /xq, one can show using a general 
solution of the system of linear differential equations, that 

B{t) = C{t)Bo + a{t)A + uj{t)Mi/2, (23) 
B'{t) = Bo + a^{t)A^ + u;^{t)M^/2^, 

where the functions ^(t), ^1(^)5 Ci't)^ <^{'t)j <^i{'t)j '^l^), ^-iid "-^il^), which determine the 
evolution of the fundamental parameters, mainly depend on a choice of Yukawa constants 
on the grand unification scale and do not depend on the initial values of the soft SUSY 
breaking parameters, /iq, and /Iq. 
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Table 1 . Mass spectrum of superpartners of observable particles and Higgs bosons for 
A^(0) = h^{0) — 10 and HeS > depending upon the choice of fundamental parameters 
A, mo, and M1/2 (all parameters and masses are given in GeV). 
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Table 2. Mass spectrum of superpartners of observable particles and Higgs bosons for 
A^(0) = h^{0) — 10 and HeS < depending upon the choice of fundamental parameters 
A, mo, and M1/2 (all parameters and masses are given in GeV). 
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Table 3. Mass spectrum of superpartners of observable particles and Higgs bosons for 
A — mo — 0, but for different initial values /i^(0),A^(0) 
(all parameters and masses are given in GeV). 
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Figure captions 



Fig.l. The lightest Higgs boson mass in modified NMSSM as a function of z — ji' /I TeV 
for hl{Q) = A2(0) = 10, mo = 0, Afg = 1 TeV, and for fx^s < 0. Thick and thin curves 
meet calculations in one-loop and two-loop approximations, respectively. 

Fig. 2. Mass spectrum in modified NMSSM as a function of z = fi'/l TeV and 
X = A/ Ml/2 for /i2(0) = A2(0) = 10, mo = 0, M3 = 1 TeV, and for /^eff < 0. Thick and 
thin curves give the masses of heavy CP-even Higgs bosons. Dotted and dashed curves 
give the masses of CP-odd Higgs bosons. The dashdot curve gives the mass of the 
heaviest neutralino. 

Fig. 3. The lightest Higgs boson mass in modified NMSSM as a function of z = n'/l TeV 
for /i2(0) = A2(0) = 10, mo = 0, M3 = 1 TcV, and for /ies > 0. Thick and thin curves 
meet calculations in one-loop and two-loop approximations, respectively. 

Fig. 4. Mass spectrum in modified NMSSM as a function of z = fi'/l TeV and 
X = A/Mi/2 for /i2(o) ^ A2(0) = 10, mo = 0, M3 = 1 TeV, and for /^eff > 0. Thick and 
thin curves give the masses of heavy CP-even Higgs bosons. Dotted and dashed curves 
give the masses of CP-odd Higgs bosons. The dashdot curve gives the mass of the 
heaviest neutralino. 

Fig. 5. Upper bound on the mass of the lightest Higgs boson in the MSSM (thin curve) 
and in the modified NMSSM (thick curve) as a function of tan (3 for M3 = 2 TeV. 
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